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DESCRIPTION 

TOMOGRAPH 
Technical Field 

5 The present invention relates to a tomograph, or more 

particularly, to a technology for reducing artifacts in a 
tomographic image attributable to a change in radiation quality 
occurring in a subject or the non-linearity in the input-output 
characteristic of a radiation detector. 

10 

Background Art 

In tomographs, a reconstructed tomographic image of a 
subject (hereinafter a reconstructed image) may contain 
artifacts. The artifacts include a ring artifact that is an 

15 annular pattern contained in a reconstructed image, and a 
dark-band artifact that is a black band-like pattern appearing 
between images of high-radiation absorbent substances. The 
major factors of the ring artifact are that the sensitivity 
of a radiation detector varies from pixel to pixel or the 

20 input-output characteristic thereof is not ideal but is 
non-linear. Moreover, the major factor of the dark-band 
artifact is that radiation quality changes in a subject (beam 
hardening) . Various methods have been proposed in efforts to 
minimize the artifacts. 

25 The most typical method of minimizing artifacts is an air 



calibration method. According to the air calibration method, 
an air image produced without a subject is produced. A 
distribution of signal intensities represented by an air image 
is equivalent to a product of a distribution of intensities 
of radiation incident on a detector by a distribution of 
sensitivities exhibited by the detector. Consequently, a 
projection image signal detected at each of pixels or pixel 
locations in the detector is divided by an air image signal, 
whereby the distribution of radiation intensities or the 
variance of detector sensitivities can be corrected. 

As another example of the artifact minimization method, 
a water correction method has been proposed (refer to, for 
example, Japanese Patent Application Laid-Open No. 7-171145) . 
The water correction method is a method classified into an 
extension of the air calibration method, and employs a water 
image, which is produced by scanning a water phantom shaped 
like a cylinder or elliptic cylinder, on behalf of the air image. 
In a case where the radiation absorption dose of the water phantom 
is close to that of a subject, the water correction method not 
only provides the same effect of correction as that of the air 
calibration method but also removes the ring artifact 
attributable to the non-linear input-output characteristic of 
a detector or the dark-band artifact attributable to the beam 
hardening. 

As still another example of the artifact minimization 



method, a phantom calibration method has been proposed (refer 
to, for example, Japanese Patent Publication No. 61-54412). 
The phantom calibration method is a method of correcting 
projection data items of a subject using a predefined transform 
function. The transform function is a polynomial expression 
for transforming measured values, which are represented by 
projection images, into theoretical values. The relationship 
between the measured values and theoretical values is drawn 
out in advance using calibration phantoms. The calibration 
phantoms having various diameters are scanned in order to detect 
the relationship between the measured values and theoretical 
values. Consequently, signal intensities represented by the 
projection images can be calibrated over a wide dynamic range 
offered by a detector. 

Disclosure of the -Invention. 

The air calibration method can correct a distribution of 
intensities of radiation incident on a radiation detector or 
a variance of sensitivities of the radiation detector, and can 
greatly minimize the ring artifact. However, the air 
calibration method cannot correct the non-linear input-output 
characteristic of the radiation detector and therefore cannot 
completely remove the ring artifact. Moreover, the dark-band 
artifact attributable to beam hardening cannot be removed. 

The water correction method not only provides the same 
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effect of correction as that of the air calibration method but 
also has the merit of removing the ring artifact and dark-band 
artifact attributable to the non-linear input-output 
characteristic of a detector. However, since the radiation 
5 absorption dose of the water phantom cannot be exactly agreed 
with that of a subject, the precision in removal is low. 

The phantom calibration method can calibrate signal 
intensities, which are represented by projection images, over 
a wide dynamic range offered by a detector. Compared with the 

10 water correction method, the phantom calibration method has 
the merit of highly precisely removing the ring artifact and 
dark-band artifact attributable to the non-linear input-output 
characteristic of a detector. However, the number of samples 
of measured values must be increased in order to improve the 

15 precision in approximation of a polynomial expression. In 
reality, it is nearly impossible to increase the number of 
samples. Only four or five samples are used in practice. The 
precision in approximation of the polynomial expression is 
therefore low. A reconstructed image having undergone the 

20 phantom calibration contains remaining artifacts. 

The reason why it is nearly impossible to increase the 
number of samples is that it is time-consuming and 
labor-intensive to acquire measured data items. For example, 
assuming that a head calibration phantom (shaped to be circular 

25 in a transaxial direction) is scanned as it is in the related 



art (Patent Document 2), only one sample of measured data is 
acquired using one calibration phantom having a certain size. 
Assuming that a chest calibration phantom (shaped to be 
non-circular in the transaxial direction) is scanned as it is 
5 in the related art, two samples of measured data items are 
acquired using one calibration phantom having a certain size 
(by scanning the chest calibration phantom in a major-axis 
direction and a minor-axis direction) . Consequently, in order 
to acquire many samples, many calibration phantoms having 

10 different sizes must be scanned in order to acquire measured 
data items. It is time-consuming and labor-intensive. 

Moreover, in order to acquire the measured data items, 
the center position of a calibration phantom must be precisely 
aligned with the center of rotation of a scanning system. It 

15 is time-consuming to position the calibration phantom. 

Furthermore, according to the conventional phantom calibration 
method, the number of samples detected at both ends of a detector 
is smaller. 

An object of the present invention is to provide a 
20 tomographic technology for improving the quality of a 

reconstructed image, which is produced by a tomograph, while 
highly precisely minimizing artifacts appearing . in the 
reconstructed image . 

The object of the present invention and the novel features 
25 thereof will be apparent from the description of the 
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specification and the appended drawings. 

In order to accomplish the above object, a tomograph in 
accordance with the present invention has features described 
below. Typical examples of the present invention will be 
described below. 

(1) A tomograph includes: a scanning system including a 
generation means for generating radiation to be irradiated to 
a subject, anda detectionmeans that is opposedto the generation 
means in order to detect the radiation transmitted by the 
subject; and a rotation means for rotating the scanning system 
about the subject. The tomograph reconstructs a tomographic 
image of the subject on the basis of a plurality of transmitted 
images produced by projecting the radiation from a plurality 
of rotational angular positions while rotating the scanning 
system about the subject. Herein, the tomograph further 
includes: a first storage means in which measured images that 
are three or more transmitted images produced by rotating the 
scanning system about a plurality of phantoms including at least 
one phantom whose section perpendicular to an axis of rotation 
of the scanning systemhas different dimensions in two directions 
orthogonal to the axis of rotation are stored; a production 
means for producing calculated images as the transmitted images 
through calculation; a second storage means in which the produced 
calculated images are stored; and a correction means for 
correcting intensities, which are represented by the 
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transmittance images of the subject, according to the measured 
images and calculated images. 

(2) In the tomograph set forth in item (1), the section 
of the at least one phantom perpendicular to the axis of rotation 

5 is substantially elliptic. 

(3) In the tomograph set forth in item (1) or (2), the 
section of at least one of the plurality of phantoms 
perpendicular to the axis of rotation is shaped substantially 
like a circle whose center is not aligned with the axis of 

10 rotation. 

(4) A tomograph includes: a scanning system including a 
generation means for generating radiation to be irradiated to 
a sub j ect and a detection means that is opposed to the generation 
means in order to detect the radiation transmitted by the 

15 subject; and a rotation means for rotating the scanning system 
about the subject. The tomograph reconstructs a tomographic 
image of the subject using a plurality of transmitted images 
produced by projecting the radiation from a plurality of 
rotational angular positions while rotating the scanning system 

20 about the subject. Herein, the tomograph further includes: 
a first storage means in which measured images that are three 
or more transmitted images produced by rotating the scanning 
system about a plurality of phantoms including at least one 
phantom whose section perpendicular to an axis of rotation of 

25 the scanning system is shaped substantially like a circle whose 
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center is not aligned with the axis of rotation are stored; 
a production means for producing calculated images as the 
transmitted images through calculation; a second storage means 
in which the produced calculated images are stored; and a 
5 correction means for correcting intensities , which are 

represented by the transmitted images of the subject, according 
to the measured images and calculated images. 

(5) In the tomograph set forth in item (4), the section 
of at least one of the plurality of phantoms perpendicular to 

10 the axis of rotation is substantially elliptic. 

(6) In the tomograph set forth in item (4), the section 
of at least one of the plurality of phantoms perpendicular to 
the axis of rotation is shaped substantially like a circle whose 
center is substantially aligned with the axis of rotation. 

15 (7) The tomograph set forth in item (1) or (4) further 

includes a phantom position calculation means for calculating 
the center position of a section of a phantom and an inclination 
of the phantom with respect to a direction parallel to the section 
according to the tomographic image of the phantom reconstructed 

20 based on the measured images. The production means determines 
a direction of projection, in which the radiation is projected 
in order to produce the calculated images, according to the 
center position and inclination. 

(8) In the tomograph set forth in item (7), the phantom 

25 position calculation means calculates the center position of 
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a section of a phantom according to the barycentric position 
in a distribution of signal intensities represented by a 
tomographic image of the phantom. 

(9) In the tomograph set forth in item (7), the phantom 
5 position calculation means calculates the inclination of a 
phantom with respect to a direction parallel to a section of 
the phantomby performing linear approximation on a distribution 
of signal intensities represented by a tomographic image of 
the phantom. 

10 (10) A tomograph includes: a scanning system including 

a generation means for generating radiation to be irradiated 
to a subject and a detection means that is opposed to the 
generation means in order to detect the radiation transmitted 
by the subject; and a rotation means for rotating the scanning 

15 system about the subject. The tomograph reconstructs a 

tomographic image of the sub j ect using a plurality of transmitted 
images produced by projecting the radiation from a plurality 
of rotational angular positions while rotating the scanning 
system about the subject. Herein, the tomograph further 

20 includes: a first storage means in which measured images that 
are three or more transmitted images produced by rotating the 
scanning system about a plurality of phantoms including at least 
one phantom whose section perpendicular to an axis of rotation 
of the scanning systemhas different dimensions in two directions 

25 orthogonal to the axis of rotation are stored; a production 



means for producing calculated images as the transmitted images 
through calculation; a second storagemeans in which the produced 
calculated images are stored; a parameter value derivation means 
for deriving parameter values to be assigned to an approximation 
function by fitting the approximation function close to the 
relationship between the signal intensities represented by the 
measured images and the signal intensities represented by the 
calculated images; a third storage means in which the parameter 
values are stored; and a correction means for correcting the 
intensities, which are represented by the transmitted images 
of the subject, according to the measured images and calculated 
images . 

(11) A tomograph includes: a scanning system including 
a generation means for generating radiation to be irradiated 
to a subject and a detection means that is opposed to the 
generation means in order to detect the radiation transmitted 
by the subject; and a rotation means for rotating the scanning 
system about the subject. The tomograph reconstructs a 
tomographic image of the sub j ect using a plurality of transmitted 
images produced by projecting the radiation from a plurality 
of rotational angular positions while rotating the scanning 
system about the subject. Herein, the tomograph further 
includes: a first storage means in which measured images that 
are three or more transmitted images produced by rotating the 
scanning system about a plurality of phantoms including at least 
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one phantom whose section perpendicular to an axis of rotation 
of the scanning system is shaped substantially like a circle 
whose center is not aligned with the axis of rotation are stored; 
a production means for producing calculated images as the 
5 transmittance images through calculation; a second storage 
means in which the produced calculated images are stored; a 
parameter value derivation means for deriving parameter values 
to be assigned to an approximation function by fitting the 
approximation function close to the relationship between the 

10 signal intensities represented by the measured images and the 
signal intensities representedby the calculated images ; a third 
storage means in which the parameter values are stored; and 
a correction means for correcting the intensities, which are 
represented by the transmitted images of the subject, according 

15 to the measured images and calculated images. 

The present invention has the advantages described below. 

(a) The number of samples of measured data items that is 
four or five according to the conventional phantom calibration 
method is increased to range from six to several thousands. 

20 Therefore, the precision in correction to be attained by the 
phantom calibration method can be improved. This leads to 
improvement of the quality of a reconstructed image. 

(b) Since the radiation attenuation constant of a phantom 
can be dynamically varied in a rotation direction in which the 

25 scanning system is rotated, measured data items falling within 
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a wide dynamic range can be acquired. 

(c) Since the radiation attenuation constant of a typical 
phantom whose section is circular can be dynamically varied 
in the rotation direction, measured data items falling within 

5 a wide dynamic range can be acquired. 

(d) Since a phantom need not be precisely positioned 
relative to the scanning system, the phantom can be readily 
placed. The work efficiency in proceeding with acquisition 
of measured data items improves. 

10 (e) Since the center position of a cross section of a phantom 

can be readily precisely detected, the precision in calibration 
to be attained by the phantom calibration method improves. 

(f ) Even when the cross section of a phantom is not circular, 
the inclination of the phantom with respect to a direction 

15 parallel to the section can be readily precisely detected. 
Therefore, the precision in calibration to be attained by the 
phantom calibration method improves. 

(g) Signal intensities represented by any transmitted 
images produced during scanning of a subject can be corrected 

20 according to the phantom calibration method. 

Brief Description of the Drawings 

Fig . 1 is a front view illustrating a tomograph in accordance 
with the first embodiment of the present invention; 
25 Fig. 2 is an explanatory diagram showing a preprocessing 
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means included in the tomograph in accordance with the first 
embodiment of the present invention; 

Figs . 3 include explanatory diagrams showing the structures 
of data items stored in the preprocessing means included in 
5 the tomograph in accordance with the first embodiment of the 
present invention; 

Fig. 4 is an explanatory diagram showing a correction table 
creating means included in the tomograph in accordance with 
the first embodiment of the present invention; 
10 Fig. 5 is an explanatory diagram concerning a correction 

table creationmethod implemented in the tomograph in accordance 
with the first embodiment of the present invention; 

Fig. 6 is an explanatory diagram showing a signal intensity 
correction means included in the tomograph in accordance with 
15 the first embodiment of the present invention; 

Figs. 7 include explanatory diagrams concerning a 
calculation method implemented in a simulation image generating 
means included in the tomograph in accordance with the first 
embodiment of the present invention; 
20 Fig. 8 is an explanatory diagram showing a calibration 

phantom position detecting means included in the tomograph in 
accordance with the first embodiment of the present invention; 

Fig. 9 is an explanatory diagram concerning a method of 
positioning a plurality of elliptic phantoms in the tomograph 
25 in accordance with the first embodiment of the present invention; 
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Figs. 10 include explanatory diagrams concerning a 
calculation method implemented in a simulation image generating 
means included in a tomograph in accordance with the second 
embodiment of the present invention; 

Fig. 11 is an explanatory diagram concerning a method of 
positioning a plurality of cylindrical phantoms in the tomograph 
in accordance with the second embodiment of the present 
invention; 

Figs. 12 include explanatory diagrams showing an example 
of another calibration phantom to be disposed in the tomograph 
in accordance with the embodiment of the present invention; 
and 

Figs. 13 include explanatory diagrams showing an example 
of an effect of image improvement attained in the tomograph 
in accordance with the first embodiment of the present invention . 

Best Mode for Carrying out the Invention 

Embodiments of the present invention will be described 
below in conjunction with the drawings. 
20 (First Embodiment) 

Fig . 1 is a front view illustrating a tomograph in accordance 
with the first embodiment of the present invention. The 
tomograph in accordance with the first embodiment includes an 
X-ray tube 1, an X-ray detector 2, a rotary panel 4, a driving 
25 motor 5, a driving belt 6, a gantry 7, a scan control means 
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100, a preprocessing means 111, a measured image memory 101, 
a correction table creating means 102, a simulation image memory 
103, a signal intensity correction means 104, an image 
reconstruction means 107, a calibration phantom position 
5 detecting means 108, an image display means 109, and a console 
110. 

Hereinafter, the assembly of the X-ray tube 1 and X-ray 
detector 2 shall be called a scanning system. The scanning 
system is fixed to the rotary panel 4. The driving motor 5 

10 rotates the rotary panel 4 and entire scanning system via the 
driving belt 6. The scanning system irradiates X-rays to a 
subject 3 from all around the subject, and produces X-ray 
transmitted images of the subject. An axis of rotation of the 
rotary panel 4 (not shown) shall be regarded as a Z axis. 

15 Moreover, coordinate axes extending horizontally or vertically 
from an origin that is the center of rotation O of the rotary 
panel 4 shall be regarded as X and Y axes respectively J An 
XYZ coordinate system defined with the X, Y, and Z axes is a 
orthogonal coordinate system. 

20 Referring to Fig. 1, a distance of an X-ray focal spot 

S in the X-ray tube 1 from the center of rotation O is typically 
690 mm. Moreover, a distance of an X-ray incidence surface 
of the X-ray detector 2 from the center of rotation O is typically 
380 mm. A time required for the rotary panel 4 to make one 

25 rotation is typically 0.5 sec. 
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A known multi-slice X-ray detector formed with ceramic 
scintillators and photodiodes is adopted as the X-ray detector 
2. The X-ray detector 2 is comprised of numerous detector 
elements (not shown) . The number of detector elements 
5 juxtaposed in the XY plane direction (a direction of channels) 
is 896, and the number of detector elements juxtaposed in the 
Z-axis direction (a direction of slices) is 32. The detector 
elements are disposed in the form of an arc substantially 
equidistantly from the X-ray focal spot S. The incidence 

10 surface of each detector element is typically 1 mm long in both 
the directions of channels and slices. The number of 
projections to be achieved during one rotation of the scanning 
systemis typically 900 . Every time the rotary panel 4 is rotated 
0.4°, one projection is performed. 

15 Next, actions to be performed in the tomograph in accordance 

with the first embodiment will be described below. The 
tomograph supports two imaging modes, that is, an actual imaging 
mode and a calibration imaging mode. An operator uses the 
console 110 to select either of the actual imaging mode and 

20 calibration imaging mode. In Fig. 1, arrows drawn with dashed 
lines indicate a flow of data treated in the actual imaging 
mode, while arrows drawn with solid lines indicate a flow of 
data treated in the calibration imaging mode. 

The actions to be performed in the tomograph in the actual 

25 imaging mode will be described below. First, an operator uses 
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the console 110 to instruct initiation of a scan. The scan 
control means 100 initiates the rotation of the rotary panel 
4 via the driving motor 5. When the rotation of the rotary 
panel 4 enters a constant-speed state, the scan control means 
100 designates the timing of irradiating X-rays from the X-ray 
tube 1 and the timing of projecting transmitted X-rays on the 
X-ray detector 2, and acquires projection data items from all 
around the subject 3. Thereafter, the preprocessing means 111 
performs preprocessing, which includes offset correction, air 
calibration, and logarithmic conversion, on the pro j ection data 
items according to a method that will be described later. The 
projection data items having undergone the preprocessing 
(hereinafter, measured images) are stored in the measured image 
memory 101. Thereafter, the signal intensity correction means 
104 reads the measured images from the measured image memory 
101, and corrects the signal intensities represented by the 
measured images according to a method that will be described 
later. The signal intensity correction means 104 references 
a correction table stored in advance in the correction table 
memory 105 so as to perform predetermined correction using the 
correction table . The correction table will be detailed later. 
Thereafter, the image reconstruction means 107 uses a known 
technology to reconstruct a tomographic image of the subject 
3 according to the output values sent from the signal intensity 
correction means 104. The reconstructed tomographic image is 
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finally displayed on the image display means 109. 

Next, actions to be performed in the tomograph in the 
calibration imaging mode will be described . In the calibration 
imaging mode, a calibration phantom to be described later is 
5 disposed on behalf of the subject 3. First, the calibration 
phantom is scanned according to the same procedure as the 
procedure in the actual imaging mode, and projection data items 
are preprocessed. The preprocessed projection data items 
(hereinafter, calibration measured images) are stored in the 

10 measured image memory 101. Thereafter, the calibration 

measured images are read by the signal intensity correction 
means 104. The signal intensity correction means 104 checks 
the correction table memory 105 to see if the correction table 
memory 105 holds a correction table. If the correction table 

15 is present, the correction table is referenced in order to 
correct the signal intensities represented by the calibration 
measured images. If the correction table is absent, the 
correction of the signal intensities represented by the 
calibration measured images is omitted. Thereafter, the image 

20 reconstruction means 107 uses a known technology to reconstruct 
a tomographic image of the calibration phantom according to 
the output values provided by the signal intensity correction 
means 104. Thereafter, the calibration phantom position 
detecting means 108 uses the tomographic image to calculate 

25 the position of the calibration phantom on the XY plane according 
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to a method to be described later. Thereafter, the simulation 
image generating means 106 uses the calculated position of the 
calibration phantom to calculate theoretic values (hereinafter 
simulation images or calculated images) from the projection 
5 data items acquired from the calibration phantom, and stores 
the results of the calculation in the simulation image memory 
103. Thereafter, the correction table creating means 102 uses 
the calibration measured images stored in the measured image 
memory 101 and the simulation images stored in the simulation 

10 image memory 103 to create correction table data, based on which 
the signal intensities represented by measured images are 
converted into theoretical values, according to a method to 
be described later, and stores the results of creation in the 
correction table memory 105. If an old correction table data 

15 is stored in the correction table memory 105, the old correction 
table data is overwritten with the new correction table data. 

Fig. 2 is an explanatory diagram showing the preprocessing 
means included in the tomograph in accordance with the first 
embodiment of the present invention. Figs. 3 are explanatory 

20 diagrams showing the structure of data items to be treated by 
the preprocessing means 111 included in the tomograph in 
accordance with the first embodiment of the present invention. 
Referring to Fig. 2 and Figs. 3, a procedure to be followed 
by the preprocessing means 111 will be described below. 

25 The X-ray detector 2 is a multi-slice detector, and has 
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896 detector elements juxtaposed in the direction of channels 
and 32 detector elements juxtaposed in the direction of slices. 
Hereinafter, N and M shall denote the numbers of detector 
elements juxtaposed in the respective directions of channels 
5 and slices. As mentioned above, 900 projections are taken 
during one rotation of the scanning system. Hereinafter, K 
shall denote the number of pro jections . Moreover, Inm(k) shall 
denote a detection signal produced by a detector element 
belonging to the n-th channel and the m-th slice (n ranges from 

10 1 to N, and m ranges from 1 to M) in the X-ray detector 2 during 
the k-th projection (k ranges from 1 to K) . 

The preprocessing means 11 performs three pieces of 
processing, that is, offset image creation, air image creation, 
and air calibration. The offset image creation and air image 

15 creation are designated for an offset image scan or air image 
scan that precedes a scan of the subject 3. The air calibration 
is designated for the scan of the subject 3. The three pieces 
of processing will be described below. 

The offset image creation is to create an average image 

20 using K offset images produced by performing an offset image 
scan (a scan not accompanied by irradiation from the X-ray tube 
1). Every time one scan is performed, the offset images are 
sequentially written over old data items in the frame memory 
200. The frame memory 200 has a data structure shown in Fig. 

25 3A, and has projection data items that number a product of N 



by M and that are equivalent to one frame data to be provided 
by the X-ray detector 2 . As soon as the offset images are stored 
in the frame memory 200, an arithmetic averaging means 201 reads 
the offset images one after another and averages th«em according 
to the formula (1) . 

A. £_] 

. - . formula ( 1 ) 

The average offset image produced by the arithmetic 
averaging means 201 is stored in the offset image memory 204. 
The offset image memory 204 has a data structure shown in Fig. 
3B, and has average offset images that number a product of N 
by M and that are equivalent to one frame data to be provided 
by the X-ray detector 2 . 

The air image creation is to create an average image using 
K air images produced by performing an air image scan (a scan 
performed without the subject 3 by irradiating X-rays from the 
X-ray tube 1) . The air images are sequentially written over 
old images in the frame memory 200. As soon as the air images 
are stored in the frame memory 200, the arithmetic averaging 
means 201 reads the air images sequentially and averages them 
according to the formula (2) . 

<C=-^IX(*) {n = \~-N,m = \~M) 

k=1 ... formula (2) 

As soon as the arithmetic averaging means 201 calculates 
an average air image, the offset correction means 202 reads 
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the average air image and compensates an offset according to 
the formula (3) . 

. . . formula ( 3 ) 
In the above arithmetic operation, an average offset image 
5 stored in the offset imagememory 204 is referenced. The average 
air image having an offset compensated by the offset correction 
means 202 is stored in the air image memory 205. The air image 
memory 205 has a data structure shown in Fig. 3C, and has average 
air images that number a product of N by M and that are equivalent 
10 to one frame data to be provided by the X-ray detector 2. 

The air calibration is performed on projection images of 
the subject in order to correct the spatial distribution of 
X-ray energy spectrum radiated from the X-ray tube 1 or the 
distribution of sensitivities exhibited by the X-ray detector 
15 2. Every time one scan is performed the projection images are 
successively written on old images in the frame memory 200. 
As soon as the projection images are stored in the frame memory 
200, the offset correction means 200 reads the projection images 
from the frame memory 200 and compensates an offset of each 
20 image compensated according to the formula (4) . 

. . . formula (4) 

The offset correction means 202 reads the projection data 
items directly from the frame memory 200 but does not read the 
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projection data items via the arithmetic averaging means 201. 
In the above arithmetic operation, the average offset images 
stored in the offset image memory 204 are referenced. As soon 
as the offsets of the projection images are compensated by the 
offset correction means 202, the air calibration means 203 reads 
the projection images having the offsets compensated, and 
performs air calibration according to the formula (5) . 

. . . formula (5) 
In the above arithmetic operation, the average air images 
stored in the air image memory 205 are referenced. The 
projection images having been calibrated by the air calibration 
means 203 are stored in the measured image memory 101. The 
foregoing series of arithmetic operations performed during air 
calibration is performed every time the projection images are 
stored in the frame memory 200, and is therefore repeated K 
times . The measured image memory 101 has a data structure shown 
in Fig. 3D, and holds projection images which number a product 
of N by M by K and which are equivalent to K frame data items 
to be provided by the X-ray detector 2. 

Fig. 4 is an explanatory diagram showing the correction 
table creating means 102 included in the tomograph in accordance 
with the first embodiment of the present invention. Fig. 5 
is an explanatory diagram concerning a correction table creation 
method implemented in the tomograph in accordance with the first 
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embodiment of the present invention. Referring to Fig. 4 and 
Fig. 5, a procedure to be followed by the correction table 
creating means 102 will be described below. 

As mentioned above, the. correction table creating means 
5 102 is used in the calibration imaging mode . In the calibration 
imaging mode, a calibration phantom to be described later is 
disposed on behalf of the subject 3. Projection data items 
of the calibration phantom are acquired. After the 
preprocessing means 111 performs air calibration on the 

10 projection data items, the projection data items are stored 
in the measured image memory 101. The data items stored in 
the measured image memory 101 and simulation image memory 103 
share the same data structure. All the projection data items 
(measured images ) , Jnm(k) (where n ranges from 1 to N, m ranges 

15 from 1 to M, and k ranges from 1 to K), having undergone the 
air calibration are stored in the measured image memory 101. 
On the other hand, theoretical values (simulation images or 
calculated images) J'nm(k) (where n ranges from 1 to N, m ranges 
from 1 to M, and k ranges from 1 to K) which the simulation 

20 image generating means 106 has calculated from the projection 
data items are stored in the simulation image memory . The method 
of calculating the simulation images J'nm(k) will be described 
later. The measured images and simulation images are ideally 
consistent with one another but actually inconsistent with one 

25 another. This is because radiation quality changes ina subject 
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(beam hardening) due to the non-linear input-output 
characteristic of the X-ray detector 2 and the energy spectrum 
of X-ray radiated from the X-ray tube 1. The non-linear 
relationship between the measured images and theoretical values 
5 may cause the ring artifact and dark-band artifact to appear 
in a reconstructed tomographic image and therefore must be 
corrected to be linear. The correction table creating means 
102 creates a correction table based on which the non-linear 
relationship is corrected. 

10 The relationship between the measured images Jnm(k) and 

the simulation images J' nm ( k) can be approximated by a polynomial 
expression or any other function. Fig. 5 graphically shows 
polynomial approximation. For the polynomial approximation, 
the axis of abscissas reads the values represented by the 

15 measured images J™ ( k) and the axis of ordinates reads the values 
representedby the simulation images J' nm(k) . Acurve isplotted 
based on points of all k values ranging from 1 to K. Thereafter, 
the curve is approximated by the polynomial function presented 
as the formula (6) below. 

20 

. . . formula (6) 

Incidentally, a known technology such as the least squares 
method is employed for the approximation. A predetermined 
value is adopted as the degree L of the polynomial function. 
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However, when the degree L is 1, since the polynomial function 
is a linear function, the polynomial function cannot express 
the non-linear relationship between the measured images and 
theoretical values. Consequently, the degree L must be 2 or 
5 larger. For precise approximation of the non-linear 

relationship, the degree L should preferably be 3 or larger. 
When the degree L is 3, the polynomial expression includes three 
coefficients anm(l) / anm(2) , and anm(3) according to the formula 
(6) . In order to calculate the three coefficients using the 

10 least squares method, the number of projections K must be 3 
or more. In general, the number of projections K should be 
equal to or larger than the degree L of the polynomial function. 
The coefficients anmCL), anmfL-l), etc., and anm(l) are derived 
from the foregoing processing, and listed in the form of a 

15 correction table. 

In order to create the correction table, the correction 
table creating means 102 first reads the measured images Jnm ( k) , 
which are acquired during all of the k-th projections ranging 
from the first projection to the K-th projection, and the 

20 simulation images J'nm(k) , which are produced from the measured 
images, from the measured image memory 101 and simulation image 
memory 103 respectively, and writes them in a buffer memory 
400 . Thereafter, a least squares approximation means 401 reads 
the data items from the buffer memory 400. and performs the 
.25 polynomial approximation according to the formula (6). The 
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obtained coefficients a nm (L), a nm (L"l) / etc., and a^l) are 
stored in the correction table memory 105. This sequence is 
repeated in order to treat all data items produced by all detector 
pixel locations (n,m) (where n ranges from 1 to N and m ranges 
5 from 1 to M). The least squares approximation means 401 is 
realizedby software installed in a dedicated or general-purpose 
computing device. 

Fig. 6 is an explanatory diagram concerning the signal 
intensity correction means 104 included in the tomograph in 

10 accordance with the first embodiment of the present invention. 
The preprocessing means 111 performs air calibration on 
pro j ection images of the sub j ect 3 acquired in the actual imaging 
mode. The resultant projection images are stored in the 
measured image memory 101. The air calibration is performed 

15 every time projection images are detected by the X-ray detector 
2. The projection data items J nm (k) (where n ranges from 1 to 
. N and m ranges from 1 to M) produced by performing air calibration 
on projection images equivalent to the k-th frame (where k ranges 
from 1 to K) are written in the measured image memory 101 . Every 

20 time the writing is completed, a polynomial calculation means 
601 reads the coefficients a nm (L) f a nm (L-l), etc., and anmfl) 
of a polynomial expression, which are associated with data items 
detected at each detector pixel location (n,m), from the 
correction table memory 105, and stores them in a buffer memory 

25 602. Moreover, the polynomial calculation means 601 reads the 
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projection data items Jnm(k) that have undergone the air 
calibration, and the coefficients a nni (L) / anm(L-l), etc., and 
anm(l) of a polynomial expression from the measured image memory 
101 and buffer memory 602 respectively, and solves the formula 
5 (6) to work out fnm(Jnm(k) ) for the purpose of signal intensity 
correction. The signal intensity correction is repeatedly 
performed in order to correct signal intensities represented 
by projection data items detectedby all detector pixel locations 
(n,m) (where n ranges from 1 to N and m ranges from 1 to M) . 
10 The results are transferred to an image reconstruction means 
107. The polynomial calculation means 601 is realized by 
software installed in a dedicated or general-purpose computing 
device . 

Figs. 7 are explanatory diagrams concerning a calculation 
15 method implemented in the simulation image generating means 
106 included in the tomograph in accordance with the first 
embodiment of the present invention. In particular, a method 
of producing simulation images in a case where an elliptic 
phantom 7 00 is adopted as the calibration phantom will be 
20 described in conjunction with Fig. 7A and Fig. 7B. An orthogonal 
coordinate system XYZ shown in Fig . 7A is a stationary coordinate 
system defined in the gantry 7 . The X-ray focal spot S is rotated 
on the XY plane, and the center of rotation is consistent with 
the origin O of the XYZ coordinate system. The XY plane meets 
25 the X-ray detector 2 along a line of intersection 702. The 
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elliptic phantom 700 is shaped like an elliptic cylinder and 
disposed so that the longitudinal direction thereof will 
substantially be aligned with the Z axis. The material of the 
elliptic phantom 700 is typically polyethylene but may be any 
5 other material such as an acrylic acid resin. Hereinafter, 
2a, 2b, and H shall denote the outer dimensions of the elliptic 
phantom 700 in the ma jor-axis direction, minor-axis direction, 
and longitudinal direction. The 2a, 2b, and H values are 
typically 350 mm, 200 mm, and 300 mm respectively. A line of 

10 intersection 701 along which the elliptic phantom 700 meets 
the XY plane depicts an ellipse. 

•As shown in Fig. 7B, a pq coordinate system is a coordinate 
system defined in the elliptic phantom 700, and the origin of 
the coordinate system is regarded as the center of an ellipse 

15 O' expressed with the line of intersection 701. Moreover, p 
and q axes shall extend in the directions of the major and minor 
axes of the substantial ellipse depicted by the line of 
intersection 701 . The elliptic phantom 700 is disposed so that 
the center point O' of the ellipse will be located near the 

20 origin O of the XYZ coordinate system and the p axis will be 
substantially aligned with the X axis . However, it is generally 
hard to dispose the phantom precisely. The position (O' x, O' y) 
of the center O' of the ellipse on the XY plane does not fully 
correspond to a point (0,0). Moreover, an angle (j) at which 

25 the p axis meets the X axis is not exactly 0. The parameter 
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(0'x,O f y) and angle § that define the position of the elliptic 
phantom700 are autonomously detected by the calibration phantom 
position detecting means 108 according to a method to be 
described later. 
5 The simulation image generating means 106 calculates 

theoretical values J'nm(k) from projection data items Jnm(k) 
(where n ranges from 1 to N and m ranges from 1 to M) that have 
undergone air calibration and are produced from the projection 
images equivalent to the k-th frame (wherein k ranges from 1 
10 toK) . The theoretical values J' nm ( k) are provided as the formula 
(7) . 

JL (k) = M p ™» m (*) in = m = 

.... formula (7 ) 

where w nm (k) denotes a distance by which the X-ray beam 703 that 
is radiated from the X-ray focal spot S during the k-th pro j ection 

15 and falls on a detector pixel location (n,m) on the X-ray detector 
2 has passed through the elliptic phantom 700. Moreover, \i p 
denotes an X-ray absorption coefficient exhibited by the 
elliptic phantom 700. Assuming that y denotes an angle of 
radiation at which the X-ray beam 703 meets the direction of 

20 channels, h^ denotes the distance of the detector pixel location 
(n, m) from the line of intersection 702, d denotes the distance 
between the X-ray focal spot S and the center of rotation O, 
and D denotes the distance between the X-ray focal spot S. and 
the incidence surface of the X-ray detector 2, the passed 
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distance w nm (k) is provided as the formula (8).. 



b t p +a t q 



. . . formula (8 ) 

where S p and S q denote p- and q-coordinates defined in the pq 
5 coordinate system in order to represent the position of the 
X-ray focal spot S, and are expressed by the formulae (9) and 
(10) respectively. 



S =(dcos& k -O^cos^ + ^sin*^ -<9')sin{* 



formula (9) 



S q =^cos^ -O' x )sin0 + (dsind k -O f y )cos$ 

. . . formula (10) 

10 In the formula (8), t p and t q denote p-axis and q-axis 

components of a unit vector representing the direction of the 
X-ray beam 703, and are expressed by the formulae (11) and (12) 
respectively. 

*p = I ~ D t7 cos (^ + Y ~ 0 k ) 

. . . formula ( 11 ) 

15 ... formula (12) 

In the formulae (9) to (12) , 6 k denotes a rotational angle 
by which the X-ray focal spot S is angularly separated from 
the X axis during the k-th projection, and is expressed by the 
formula ( 13) . 
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e k =2*^ (k=i~K) 

... . formula (13) 
The simulation image generating means 106 uses the formulae 

(7) to (13) to calculate theoretical values J'nm(k) from 

projection images detected at all pixel locations (n,m) , where 

n ranges from 1 to N and m ranges from 1 to M, during all of 

the k-th projections where k ranges from 1 to K. The results 

of the calculation are stored in the simulation image memory 

103. The simulation image generating means 106 is realized 

by software installed in a dedicated or general-purpose 

computing device. 

Fig. 8 is an explanatory diagram showing the calibration 

phantom position detecting means 108 included in the tomograph 

in accordance with the first embodiment of the present invention . 

As already described, the elliptic phantom 700 is disposed so 

that the center position O' will be nearly consistent with the 

center O of the XY plane. However, the center position O' need 

not be highly precisely consistent with the center O, as long 

as a deviation of the center position O' from the center O falls 

within several centimeters. Since such rough precision is 

permitted, man-hours required for disposition of the 

calibration phantom are reduced. The calibration phantom 

position detecting means 108 autonomously detects a deviation 

of the position of the elliptic phantom 700 caused by the above 

disposition. The positional deviation is determined with the 
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center position 0' of the elliptic phantom 700 and an inclination 
§ at which the major axis of the elliptic phantom 700 (p axis) 
meets the X axis. The detected O' and <|> values are referenced 
by the simulation image generating means 106. 
5 A CT value thresholding means 800 reads a CT reconstructed 

image of the elliptic phantom 700 from the image reconstruction 
means 107. Hereinafter, signal values representing the CT 
reconstructed image shall be defined as signal values R(i, j) . 
i and j denote pixel locations in the X axis direction and the 

10 Y axis direction, respectively, in the reconstructed image 
(where i ranges from 1 to I and j ranges from 1 to J) . The 
CT value thresholding means 800 references a threshold R t 
recorded in advance in a threshold memory 8 03, and compares 
the threshold R t with the signal values R(i, j) . If R t <R(i, j) 

15 is established, the signal values R(i,j) are replaced with Is. 
If R(i,j)< R t is established, the signal values R(i,j) are 
replaced with 0s. Incidentally, the threshold R t is re-set to 
an intermediate value between a CT number exhibited by image 
data representing the internal region of the elliptic phantom 

20 and a CT number exhibited by image data representing the external 
region thereof. Consequently, the portion of a CT 
reconstruction image binary-coded by the CT value thresholding 
means 800 which shows the internal region of the elliptic phantom 
is represented by Is, the portion thereof which shows the 

25 external region thereof is represented by 0s . The binary-coded 
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CT reconstructed image is read by each of a barycenter 
calculation means 801 and a slant calculation means 802. The 
barycenter calculation means 801 is a means for calculating 
a barycentric position in the elliptic phantom 700. The 
barycentric position is consistent with the center position 
O' of the elliptic phantom 700, and is calculated according 
to the formula (14) below. 



0' = {0' x ,0' y ) = 



; j 



. . . formula (14 ) 
where Xi and Y ± denote X- and Y-coordinates defined in the XY 

10 coordinate system in order to represent the position of a pixel 
(i,j). The slant calculation means 802 is a means for 
calculating an inclination <() at which the major axis of the 
elliptic phantom 700 meets the X axis. The slant calculation 
means 802 performs linear approximation on the pixels (Xj_,Yi) 

15 represented by the signal values R(i, j) of Is according to the 
formula (15) . 



20 



Y =AX+B 



. . formula (15) 



At this time, the inclination <|> is calculated according 
to the formula (16) . 



<f> = tan" 1 A 



formula (16) 



35 



Incidentally, a known technology such as the least squares 
method is adopted for the linear approximation. Moreover, the 
calibration phantom position detecting means 108 is realized 
by software installed in a dedicated or general-purpose 
computing device. 

Figs. 13 are explanatory diagrams concerning an effect 
of improvement in image quality exerted by the tomograph in 
accordance with the first embodiment of the present invention. 
Above all, Fig. 13A illustrates how to dispose an assessment 
subject 1300. Fig. 13B and Fig. 13C show prof iles indicating 
CT numbers, which are represented by a reconstructed image in 
cases where calibration is not performed or performed, in 
relation to positions on the Y axis. The assessment subject 
1300 is a water phantom shaped like a cylinder of 350 mm in 
diameter. For a scan, the tube voltage of the X-ray tube 1 
is set to 120 kV, the tube current thereof is set to 200 mA, 
and the other scanning conditions are identical to those 
described in relation to the first embodiment . In the case 
where calibration is not performed, the profile 1301 expresses 
poor homogeneity. A maximum difference of a CT number from 
another is 59 HU . In contrast, in the case where calibration 
is performed, the profile 1302 expresses excellent homogeneity. 
The maximum difference of a CT number from another is 4.8 HU . 
Consequently, according to the present invention, the precision 
in deriving CT numbers represented by a reconstructed image 



36 



improves, and image quality improves. 

The tomograph in accordance with the first embodiment has 
been described so far. The present invention is not limited 
to the first embodiment but can be modified in various manners 
5 without a departure from the gist. For example, according to 
the first embodiment, the elliptic phantom 700 is limited to 
one size. Alternatively, as shown in Fig. 9, a plurality of 
elliptic phantoms 700a to 700d having different sizes may be 
used to create a correction table. However, the elliptic 

10 phantoms 700a to 700d shall be disposed so that the center 
positions thereof will be located near the center of rotation 
O in the scanning system. Moreover, the elliptic phantoms 700a 
to 700d shall be disposed so that they will fully fall within 
a field of view 900 determined by the X-ray detector 2. 

15 Furthermore, the correction table creating means 102 performs 
polynomial approximation on calibration measured images, which 
are acquired from the elliptic phantoms 700a to 700d, and 
simulation images according to the formula (6) , and stores the 
calculated coefficients anm(L), anm(L-l), etc., and anm(l) in 

20 the correction table memory 105. 
(Second Embodiment) 

A tomograph in accordance with the second embodiment of 
the present invention will be described below. The tomograph 
in accordance with the second embodiment of the present invention 
25 uses a cylindrical phantom 1000 on behalf of the elliptic phantom 
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700 employed as a calibration phantom in the first embodiment. 
The other components of the tomograph are identical to those 
of the first embodiment described in conjunction with Fig. 1 
to Fig. 6. An iterative description will be omitted. 

Figs . 10 are explanatory diagrams concerning a calculation 
method implemented in the simulation image generating means 
106 included in the tomograph in accordance with the second 
embodiment of the present invention. Referring to Figs. 10, 
a simulation image generation method will be described on the 
assumption that the cylindrical phantom 1000 is adopted as a 
calibration phantom. As already described, the orthogonal 
coordinate system XYZ is a stationary coordinate system defined 
in the gantry 7. The X-ray focal spot S is rotated on the XY 
plane, and the center of rotation is consistent with the origin 
O of the XYZ coordinate system. The XY plane meets the X-ray 
detector 2 along a line of intersection 1002. The cylindrical 
phantom 1000 has a cylindrical shape and is disposed so that 
the longitudinal direction thereof will substantially be 
aligned with the Z axis. The cylindrical phantom 1000 is made 
of a substantially homogeneous material at a substantially 
uniform density. The material of the cylindrical phantom 1000 
is typically polyethylene. Alternatively, any other material 
such as an acrylic resin may be substituted for the polyethylene . 
Hereinafter, 2r and H shall denote the outer diameter and height 
of the cylindrical phantom 1000 respectively. The 2r and H 



values are typically 250 mm and 300 mm. 

The line of intersection 1001 along which the cylindrical 
phantom 1000 meets the XY plane depicts substantially a circle, 
and the center of the circle shall be a point O' . The cylindrical 
5 phantom 1000 is disposed so that the center point O' of the 
substantial circle will.be located at a position (0'x,0'y) 
different from a point on an axis of rotation, that is, the 
origin O of the XYZ coordinate system. 

The parameter value (0'x,0'y) representing the position 

10 of the cylindrical phantom 1000 is autonomously detected by 
the calibration phantom position detecting means 108. The 
calibration phantom position detecting means 108 autonomously 
detects the center point O' according to the same method as 
the method described in con junction with Fig . 8. However, when 

15 the cylindrical phantom 1000 is an object of tomography, unlike 
when the elliptic phantom 700 is employed, the inclination (|) 
need not be detected. Consequently, calculation to be 
performed by the slant calculation means 802 is omitted. Only 
the center point O' detected by the barycenter calculation means 

20 801 is referenced by the simulation image generating means 106. 

The simulation image generating means 106 calculates 
theoretical values J'nmCk) from projection data items Jm^k) 
(where n ranges from 1 to N and m ranges from 1 to M) that are 
produced from projection images equivalent to the k-th frame 

25 (where k ranges from 1 to K) and that have undergone air 



39 



calibration. The theoretical values J'nm(k) are expressed as 
the formula (17) . 



where WnnJk) denotes a distance by which an X-ray beam 1003 that 
5 is radiated from the X-ray focal spot S during the k-th pro j ection 
and falls on a detector pixel location (n,m) on the X-ray detector 
2 passes through the cylindrical phantom 1000. Moreover, ja q 
denotes an X-ray absorption coefficient exhibited by the 
cylindrical phantom 1000. Assuming that y denotes an angle 

10 of radiation at which the X-ray beam 1003 meets the direction 
of channels, h™ denotes a distance of the detector pixel location 
(n,m) from the line of intersection 1002, d denotes a distance 
between the X-ray focal spot S and the center of rotation O, 
and D denotes a distance between the X-ray focal spot S and 

15 the incidence surface of the X-ray detector 2, the passed 
distance w nm (k) is provided as the formula (18). 



where S' x and S f y are expressed as the formulae (19) and (20) 
respectively . 

S' x =dcos0 k -O' x 



formula (17) 




. . . formula ( 18 ) 



20 



formula (19) 
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S' y =dsin0 k -O' y 

. . . formula (20) 
In the formula (18) , t x and t y denote X- and Y-axis components 
of a unit vector representing the direction of the X-ray beam 
1003, and are expressed as the formulae (21) and (22) 
5 respectively. 

- D , ^ . 

V-O- + hr nm 

. . . formula (21) 

. . . formula (22) 
In the formulae (19) to (22) , 9k denotes a rotational angle 
by which the X-ray focal spot S is angularly separated from 

10 the X axis during the k-th projection, and is expressed as the 
formula (13) . The simulation image generating means 106 uses 
the formulae (17) to (22) to calculate theoretical values J' nm ( k) 
from projection images detected at all pixel locations (n,m) , 
where n ranges from 1 to N and m ranges from 1 to M, during 

15 k-th projections where k ranges from 1 to K. The results of 
the calculation are stored in the simulation image memory 103. 
The simulation image generating means 10 6 is realized by software 
installed in a dedicated or general-purpose computing device. 
The tomograph in accordance with the second embodiment 

20 has been described so far. Needless to say, the present 

invention is not limited to the second embodiment but can be 
modified in various manners without a departure from the gist 



thereof. For example, according to the second embodiment, the 
cylindrical phantom 1000 is limited to one size . Alternatively, 
a plurality of cylindrical phantoms 1000a to lOOOd having 
different sizes as shown in Fig. 11 may be used to create a 
correction table. However, the cylindrical phantoms 1000a to 
lOOOd shall be disposed so that the peripheries thereof will 
be substantially inscribed to the field of view 900 determined 
by the X-ray detector 2 . However, all the cylindrical phantoms 
1000a to lOOOd shall be disposed to fully fall within the field 
of view 900 determined by the X-ray detector 2. Moreover, the 
correction table creating means 102 performs polynomial 
approximation on calibration measured images, which are 
acquired from all the cylindrical phantoms 1000a to lOOOd, and 
simulation images according to the formula (6), and stores 
calculated coefficients anm(L), anm(L-l), etc., and a^l) in 
the correction table memory 105. 

In the aforesaid first and second embodiments , the elliptic 
phantom 700 or cylindrical phantom 1000 is adopted as a 
calibration phantom.. However, the calibration phantom is not 
limited to these types. For example, an abdominal phantom 1200 
simulating the human abdomen as shown in Fig. 12A or a chest 
phantom 1201 simulating the human chest as shown in Fig. 12B 
may be adopted. Herein, the chest phantom 1201 has holes 1202 
and 1203 bored in order to simulate the lung fields of a human 
body. When the shape of the calibration phantom resembles the 
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shape of the subject that is an object of tomography, the number 
of scattered X-rays produced during a scan of the calibration 
phantom becomes close to the number of scattered X-rays produced 
during a scan of the subject 3. Therefore, the precision in 
5 correction to be performed by the signal intensity correction 
means 104 improves. 

In the first and second embodiments, one calibration 
phantom of a certain shape is employed. The phantoms employed 
in the first and second embodiment's respectively may be used 

10 in combination in order to acquire calibration measured images, 
and the calibration measured images and simulation images 
produced from the calibration measured images may be used to 
create the correction table 105. 

Moreover, the phantoms employed in the first and second 

15 embodiments respectively may be used in combination with various 
calibration phantoms of other different shapes or sizes 
(including a phantom whose section perpendicular to an axis 
of rotation is substantially circular and which is substantially 
aligned with the axis of rotation) in order to acquire 

20 calibration measured images, and the calibration measured 
images and simulation images produced from the calibration 
measured images may be used to create the correction table 105. 

As described so far, the present invention realizes a 
tomographic technology for acquiring many samples of measured 

25 data items from a calibration phantom by performing simple 
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measurement according to a phantom calibration method. 
Consequently, the precision in polynomial approximation based 
on the phantom calibration method improves, and the quality 
of a reconstructed image improves. 

Industrial Applicability 

According to the present invention, artifacts appearing 
in a reconstructed image produced by a tomograph are reduced 
to improve the quality of the reconstructed image. 



